Wave scattering on a domain wall in a chain of TT-symmetric couplers 
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We study wave propagation in linear arrays composed of pairs of conjugate waveguides with 
balanced gain and loss, i.e. arrays of the PT-symmetric couplers, where the linear spectrum is 
known to feature high-frequency and low- frequency branches. We introduce a domain wall by 
switching the gain and loss in a half of the array, and analyze the scattering of linear waves on 
this defect. The analysis reveals two major effects: amplification of both reflected and transmitted 
waves, and excitation of the reflected and transmitted low-frequency and high-frequency waves by 
the incident high-frequency and low-frequency waves, respectively. 



PACS numbers: 42.25.Bs, 11.30.Er, 42.82. Et, 42.81. Qb 

I. INTRODUCTION 

It is well established that non-Hermitian Hamiltonians, 
subject to the constraint of the parity-time (PT) symme- 
try, i.e., the equilibrium between spatially separated loss 
and gain, which are set as mirror images of each other, 
may give rise to entirely real spectra, provided that the 
strength of the gain and loss does not exceeds a critical 
level Although, generally speaking, PT-symmetric 
settings belong to the class of dissipative systems, they 
can support continuous families of both linear and non- 
linear modes, thus resembling the main properties of con- 
servative systems. In fact, 'PT-symmetric models lie at 
the borderline between conservative and truly dissipative 
dynamical systems. 

It is straightforward to implement PT-symmetric com- 
plex potentials, V{x), which must be subject to the afore- 
mentioned equilibrium condition, V(x) = V*(— x), by 
symmetrically juxtaposing elements accounting for the 
gain and loss [2]. This possibility was elaborated in a 
number of theoretical [l]-[8] and experimental @ stud- 
ies. 

In optics, the basic PT-symmetric element can be re- 
alized as a pair of linearly coupled waveguides, one with 
a lossy core and the other one carrying a matched com- 
pensating gain [4j. A chain composed of such coupled 
elements was introduced in Ref. [5], assuming that each 
amplified and dissipative waveguide was linearly linked 
to an adjacent waveguide of the opposite sign, belonging 
to the neighboring pair. 

Another type of such a PT-symmetric system was pro- 
posed in Ref. Q, with the gain- and loss-carrying units 
coupled to their neighboring counterparts of the same 
sign [see Fig. HJa) below] . Assuming that each unit also 
carried the conservative cubic nonlinearity, discrete soli- 
tons were found in this setting. 

The subject of this paper is a system of the general 
same type (although without nonlinearity), but includ- 
ing a defect in the form of the domain wall (DW), as 



shown below in Fig. [IJb). A natural dynamical problem 
to consider in such an array, in addition to the DW it- 
self, is the scattering of linear waves on the DW, which is 
another subject of the present work (the scattering of lin- 
ear waves on an isolated PT-symmetric complex, includ- 
ing such specific features as amplification of transmitted 
and reflected waves and Fano resonances, was analyzed 
in Refs. @ and Q). Below we demonstrate that the 
scattering of waves on the DW gives rise to nontrivial 
effects, including the transformations between different 
branches of the traveling-wave modes and their amplifi- 
cation. The fact that we consider the scattering of inci- 
dent waves with real frequencies, and the generation of 
transmitted and reflected waves which are carried by real 
frequencies too, implies that we are dealing with the case 
when the spectrum of the PT-invariant system is purely 
real, under the condition that the gain-loss coefficient is 
kept below the critical value [see Eq. below]. 

It is relevant to mention that this scattering problem 
is related to the analysis of transport and scattering pro- 
cesses governed by non-Hermitian Hamiltonians [10j . In 
the general case, such Hamiltonians are not subject to 
the PT-symmetry constraint, therefore the correspond- 
ing spectrum is complex. In particular, the analysis of 
the generic spectra demonstrates that they contain a few 
eigenvalues with especially large imaginary parts. The 
rapid decay of the corresponding eigenstates may be con- 
sidered as an analog of superradiance in optics and, 
naturally, it strongly affects dynamical features of such 
systems. It remains to understand if PT-symmetric sys- 
tems may give rise to similar "superradiant" states in the 
respective complex spectrum (above the transition from 
the real spectrum). 

The paper is organized as follows. Section [TT] intro- 
duces the model, and also it discusses the geometry of a 
domain-wall defect introduced into the chain. Analytical 
results for the transmission and reflection coefficients are 
presented in Sec. lIIIl whereas the dependence of the scat- 
tering coefficients on the system parameters is analyzed 



in Sec. IIV1 Finally, Sec. IVl concludes the paper. 



II. THE MODEL 

We consider the array of paired waveguides with bal- 
anced gain and loss (7), i.e., "PT-symmetric couplers, as 
shown in Fig. [Ha). The linear-coupling constant for the 
waveguides in the VT elements is normalized to be 1, 
while the coefficient of the linear coupling between adja- 
cent elements in the array is C±. Thus, the model is based 
on the following system of linear Schrodinger equations: 
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where u n (z) and v n (z) are the complex amplitudes in 
the amplified and damped waveguides at each site of the 
array. Actually, Eq. (JTJ) is a linearized version of the 
model introduced in recent work [8] . 

The DW in the array is created by switching the gain 
and loss in a half of the chain; generally, the constant 
of the linear coupling between the halves, C2, may be 
different from the regular value, C\ [see Fig. [If b)] . Thus, 
the array with the embedded DW is described by the 
following equations: 
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7U0 + iv + iC\(u-i - u ) + iC 2 (ui - u ), 
-jv Q + iu + id(v-i - v ) + id(vi - v ), (4) 
— 7U1 + ivi + iCi(u 2 - ui) + iC 2 {u a — ui), 
7«i + iui + iCi(v 2 - vi) + iC 2 {v — vi), (5) 

III. ANALYTICAL RESULTS 



First, we consider the wave propagation in the array 
without the DW. The corresponding solution to Eq. (TT]) 
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Figure 1: (a) A schematic of the chain of PT-symmetric cou- 
plers composed of waveguides with gain and loss, which are 
designated by dark and bright circles, respectively, (b) The 
model with the domain wall, created by inverting gain and 
loss in a half of the chain. 
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where wavenumber k may be complex, while frequency 
uj is real. Substituting Eq. ^ into Eq. ([1]), one finds the 
high-frequency (HF) and low-frequency (LF) branches of 
the dispersion relation for the linear waves, denoted by 
subscripts h and I, respectively: 
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In particular, for k = in with real n we obtain an 
exponential- wave (EW) solution to Eq. (flj, while k = 
in + n gives rise to a staggered exponential-wave (SEW) 
one. Continuous-wave (CW) solutions correspond to 
k = n, i.e., real wavenumbers. Only these three types of 
the waves admit real frequencies ui, which, in addition, 
requires that the gain/loss coefficient must be smaller 
than the constant of the coupling between the amplified 
and dissipative waveguides in each VT element, i.e., 



(9) 



An example of the spectrum determined by Eqs. ([7]) 
and ([8]) is presented in Fig. [2j In the EW region, — 00 < 
k < 0, u)h(in) and are plotted. The range of < 

« < 7r corresponds to the CW solutions, where we show 
the frequencies as uih{^) and uji(k). At k > ti, we have 
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expressions into Eqs. we obtain 
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Figure 2: Spectrum of the array without the domain wall, 
obtained from Eq. |QJ for Ci = 2, 7 = 0.5. The HF and LF 
branches, m and wj, are defined by Eqs. 0,©, with k — in 
for —00 < k < (the EW region), k — k for < k < tt (the 
CW region), and k = ir + i(n — n) for n < k < 00 (the SEW 
region) . 



the SEW region, in which the branches of the dispersion 
relation are plotted as uih {it + i [k — it] ) and lji (tt + i [n — n] ) . 

Now we proceed to the scattering of linear waves on the 
DW, within the framework of Eqs. ©-©■ To this end, 
we consider an incident wave, with frequency u>i (and the 
intensity set equal to 1), which approaches the DW from 
the left. Only the case when the incident wave belongs 
to the LF branch is studied below in an explicit form, 
as the case of the HF incident wave can be considered 
similarly. 

We look for the scattering solution to Eqs. J!])-© as 
follows: 
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for n > 1. Here the characteristic of the plane- wave 
components of the solution are sin Si — sin^ = —7, 
cos (5; = — cos Sh = yl — 7 2 , and kh, fa are to be found 
from Eqs. (J?) and Eqs. ©, respectively, where we set 
Wfc = uji = uji. 

Amplitudes i?;, 7] and i?^, T^, defined in expressions 
(1101) and are complex reflection and transmission 
coefficients for the LF and HF waves. Substituting these 
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nated in Fig. [2] Taking u;j = ujj, we conclude that both 
kh = and fc; = 



The so obtained reflection and transmission coefficients 
can be analyzed for the LF incident wave of any of the 
three types, EW, CW, or SEW. As seen from Eqs. (QUI) 
and (fTTj) . the incident wave of the LF type generates, 
generally speaking, both LF and HF reflected and trans- 
mitted wave components, whose frequencies are identical 
to cjj. 

We start by considering two examples, which are desig- 

UJl 

are in the CW region, meaning 
that the reflected and transmitted HF and LF waves are 
all of CW type (here, the asterisk does not stand for the 
complex conjugate). On the other hand, for = u>}* we 
have kh — —in h * and fa = meaning that the incident 
wave is the LF of the CW type, while the reflected and 
transmitted ones have the HF and LF components of the 
EW and CW types, respectively. 

Care should be taken as concerns the choice of the 
sign in front of imaginary part of kh and fa for the EW. 
For the plus (minus) sign, the reflected and transmitted 
waves of the EW type exponentially decrease (increase) 
with the increase of the distance from the DW. Both cases 
are physically meaningful in the VT system, but in the 
following we focus only on the evanescent (exponentially 
decaying) EWs. 

Several generic examples of the wave-intensity profiles 
are present in Fig. [3] in the array of couplers for differ- 
ent frequencies of the incident wave, wj, which are cho- 
sen with regard to the dispersion relation displayed in 
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Fig. [2j As mentioned above, the LF incident wave gen- 
erates reflected and transmitted waves belonging to the 
HF and LF branches, with the same frequency uij. The 
following four cases are present in Fig. [3] (a) uji = — 1, 
with k h = 0.931843 i, ki = 0.258102 i, hence the inci- 
dent wave is the EW of the LF type, the corresponding 
HF and LF components also being EWs [see Fig. [3] (a) , 
where the logarithmic scale is used on the vertical axis]. 
The intensity of the incident wave increases with dis- 
tance form the DW, while the reflected and transmitted 
waves are evanescent, (b) ojj = 0.5, with kh = 0.424603*, 
ki = 0.851981, hence the incident LF wave is a CW, while 
the reflected and transmitted waves include LF CW and 
HF evanescent component, the latter rapidly decaying 
with the distance from the defect [Fig. [3th)]. (c) ojj = 6, 
with kh = 1.85823, fc; = 2.36958, in which case all the 
waves are of the CW type [Fig. |3] (c)]. (d) uij = 7.5, 
with k h = 2.28961, h = n + 0.424603?, which makes the 
incident LF wave an SEW, while the reflected and trans- 
mitted waves include the HF CW and LF SEW terms 
[see Fig. [3] (d) , where the logarithmic scale is again used 
on the vertical axis] . 
Note that, at 
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the splitting between the LF and HF branches in Fig. [2] 
becomes so large that their CW regions do not have 
common frequencies. Under this condition, the case pre- 
sented in Fig.[3tc), with all the wave components being of 
the CW type, is impossible. In terms of lu 
and w max = 4Ci - y/l -7 . 
the largest frequencies of the HF and LF branches in 
the CW region (see [2]) , condition (|T6| is tantamount to 

Wmax ^ k-miin • 

In the following we take C\ > 0.5, which rules out 
condition (fT6"|) for all < 7 < 1. We thus always include 
the possibility of having all the reflected and transmitted 
waves of the CW type, excited by the incident LF CW. 

In the following Section we present the analysis of 
the transmission and reflection coefficients as functions 
of parameters of the system, assuming, as said above, 
w ln i n < w max , and focusing on the range of the incident- 
wave's frequency u> m i n < loj < cJ max , when this wave is 
of the most relevant CW type. Further, it is convenient 
to define the normalized frequency, 
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so that uj'j — at ljj — w m i n , and lj'j = 1 at u)i = oj n 



IV. TRANSMISSION AND REFLECTION 
COEFFICIENTS 

Here we analyze the transmission and reflection coef- 
ficients for the LF incident wave, given by Eq. ([12")) and 
Eq. (fT3")) , varying parameters 7, Ci, and C = Ci/CV 
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Figure 3: Examples of solutions of the scattering problem 
in the array of couplers for Ci = 2, C2 = 1, 7 = 0.5, 
and different values of the incident-wave's frequency, (a) 
The incident low-frequency exponential wave excites the high- 
frequency evanescent wave, with u>i = — 1, kh = 0.931843 i, 
ki — 0.258102 i. (b) The incident low-frequency CW excites 
a high-frequency exponential wave, with oji = 0.5, kh = 
0.424603 i, ki = 0.851981. (c) The incident low-frequency 
CW excites a high-frequency CW, with wj = 6, kh = 1.85823, 
ki = 2.36958. (d) The incident low-frequency staggered ex- 
ponential wave excites a high-frequency CW, with uji = 7.5, 
kh ~ 2.28961, ki — ir + 0.424603 i. Dots and open circles show 
\u n \ 2 and \v n \ 2 , respectively. 



As mentioned above, we consider only the case of the 
CW incident wave as the most natural one. Then, two 
possibility may be expected, with the incident LF wave 
exciting either HF-EW or HF-CW. 

It is useful to study first the case with no gain and loss, 
7 = 0. Under this condition, Eq. (TT2l) reduces to 
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From here, we immediately arrive at the conclusion that 
the LF incident wave excites only the LF reflected and 
transmitted waves, satisfying condition |i?;| 2 + |T;| 2 = 
1 due to the energy conservation. Thus, the excitation 
of the HF reflected and transmitted waves by the LF 
input is a specific effect of the VT system, due to the 
presence of the gain and loss in it. Expressions for Ri 
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Figure 4: \Ri\ 2 (the thick line) and |T ( | 2 (the thin line) as 
functions of uii in the case of the incident low-frequency CW 
for Ci = 2, C2 = 1, and 7 = (no gain and loss). In this case, 
the high-frequency waves are not excited by the low-frequency 
input (R h =Q,T h = 0). 
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and Ti in Eqs. (|18j) become singular, with a vanishing 
denominator, only when C = 1 and fc; = it, which means 
the absence of the DW in the system (in which case Ti = 1 
and Ri — 0). Figure@]displays \Ri\ 2 and |T/| 2 as functions 
of co i for parameters 7 = and C = 2. 

In the presence of the gain and loss (7 > 0) the solu- 
tion of the scattering problem for the LF incident wave 
contains Rh,Th ^ 0, i.e., the HF components are ex- 
cited. Note, however, that Rh = Th = 2/ at ki = or 
ki = it, i.e., at the borders of the considered range of the 
incident- wave's frequency, uij. Also, XJ = if k% = 
or C = 2/(1 - e~ ikh ) (the latter is possible in the SEW 
region, which is not under consideration here). 

We now fix 7 = 0.5, C\ = 2, and plot the (typical) 
results, produced by Eqs. (fl~2j) for C = 2, C = 1, and 
(7 = 0.5. " _ 

Figure [5] shows that, at C = 2, all four reflection and 
transmission coefficients have finite maxima close to the 
point lo'j = 0.001 (cjj = 0.872). Note that (a') and 
(b') show blowups of (a) and (b), respectively, in the 
vicinity of the maxima. Also, T/ vanishes at lu'j = 
(ui = 0_866025). 

For C = 1, the reflection and transmission coeffi- 
cients are presented in Fig. [5] All four coefficients di- 
verge at uj'j = 0.0022 (uj = 0.88 ) and u)\ = 0.999 
(luj = 7.1225) because at these two points the denomi- 
nator D in Eq. (|12p. which is common for all coefficients, 
vanishes. On (a'),(b') we show the blowups of (a),(b) for 
the maxima close to uj'j — 0, while in (a"),(b") for the 
maxima close to uj'j = 1. The coefficient T\ drops to zero 
at uj'j = (wj = 0.866). 

For C = 0.5, the reflection and transmission coeffi- 
cients are displayed in Fig. [7J All the four coefficients 



have finite maxima at uj'j = 0.00446 (uj = 0.894), and T; 
vanishes at uj'j = (ujj = 0.86603). 

The above results can be summarized by saying that 
for C ^ 1 all the four reflection and transmission coef- 
ficients attain finite maxima at uj'j close (but not equal) 



Figure 5: Relative intensities of the reflected and transmitted 
waves as functions of uj'j in the case of the low-frequency inci- 
dent CW, for 7 = 0.5 an C — 2. Thick and thin lines depict 
\Ri\ 2 , \R h \ 2 in (a) and \T t \ 2 , \T h \ 2 in (b), respectively, (a') 
and (b') are the blowups of (a) and (b), respectively. 
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Figure 6: Same as in Fig. [5] but for C — 1. 
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Figure 8: |7?(| 2 (thick line), |T;| 2 (thin line), \R h \ 2 (dashed 
thick line), \Th\ 2 (dashed thin line) as functions of 7 in the 
case of low-frequency incident CW for Ci = 2, C2 = 1. (a) 
ki = 2 with real fch, and (b) ki — 1 with fcj, changing from 
imaginary to real at 7 = 0.39. 



Figure 7: Same as in Figs. [5] and [6] but for C = 0.5. 



to zero. For the special case of C = 1 all four coefficients 
diverge at close to zero and close to unity. 

It is relevant to note that, if Ri = at a particular 
value of wj, then Tj = 1, and, conversely, Ri = 1 if T\ = 
0. These cases correspond to the full transmission and 
full reflection for the LF waves, respectively. However, 
at 7 ^ the HF reflection and transmission waves are 
excited too, therefore, in fact, these cases imply not full 
transmission and reflection, but rather full conversion of 
LF waves into their HF counterparts. 

It is also interesting to consider the effect of amplifi- 
cation of the reflected and transmitted waves, which is 
possible in the presence of the gain and loss, 7 > [y]. To 
this end, we fix C\ = 2 and C2 = 1 and vary the gain/loss 
parameter within the range of < 7 < 1 at fixed fcj, 
which leads to the variation of wj from 4Ci ( 1 — cos ki ) — 1 
to 4Ci(l — cosfc;). Accordingly, kh changes too, taking 
both real and imaginary values. For this case, coefficients 
|i? ; | 2 , |7]| 2 , \R h \ 2 , and \T h \ 2 are shown versus 7 in Fig.0 
for ki = 2 (a) and ki — I (b). Note that kh remains real 
in (a) for all 7, while in (b) kh changes from imaginary 
to real at 7 = 0.39. As seen in the figure, in (a) all the 
reflection and transmission coefficients increase with 7. 
On the other hand, in panel (b) they feature additional 
extrema in a vicinity of the point where kh changes from 
imaginary to real. 

The amplification of the reflected and transmitted 
waves can be defined in terms of their total intensities. 
Thus, in the case of the incident LF CW, when both re- 
flected and transmitted HF waves are of CW type, the 
amplification of reflection and transmission takes place 
when |i?/| 2 + \R h \ 2 > 1 and |T ; | 2 + \T h \ 2 > 1, respec- 
tively (as shown in Fig. |9|). We notice that the range of 
the incident-wave's frequency, ojj, where the amplifica- 
tion occurs, increases with 7. Also note that in Fig. [H] (a) 
the reflection and transmission coefficients feature sharp 
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Figure 9: |iii| 2 + \R h \ 2 and |T;| 2 + |T ;i | 2 (thick and thin lines, 
respectively) as functions of uj'j in the case of low-frequency 
incident CW for Ci = 2, C 2 = 1, and 7 = 0.01 (a), 7 = 0.5 
(b), 7 = 0.9 (c), 7 = 0.999 (d) Note that the vertical scale is 
logarithmic. 



maxima in a vicinity of tu'j — 0. In the limit 7 = 0, 
as it has been already mentioned, the system becomes 
conservative and these maxima disappear. 



V. CONCLUSIONS 

We have studied the scattering of linear waves on a 
domain wall introduced into a waveguide array com- 
posed of 'PT-symmetric waveguide pairs. Such arrays 
support the propagation of HF (high-frequency) and LF 
(low-frequency) waves. Considering incident LF waves 
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of various types (continuous waves or unstaggered and 
staggered exponential waves) , we have derived the corre- 
sponding reflection and transmission coefficients and an- 
alyzed their dependence on the system parameters. The 
case of the HF incident wave can be analyzed similarly. 

We have found that the LF incident wave generates 
both LF and HF reflected and transmitted waves, pro- 
vided that the gain and loss are present (7 > 0). We 
also demonstrated that both reflected and transmitted 
waves can be substantially amplified, provided that the 
gain 7 > is present. The range of the incident-wave's 
frequency where the amplification takes place expands 
with the increase of 7. 

Our results suggest that the use of "PT-symmetric ele- 
ments in waveguide arrays offers various possibilities for 
manipulations of optical signals in photonic lattices. It 
may be interesting to add nonlinearity to the system. In 



addition to the formation of solitons [8j, the nonlinear- 
ity may give rise to spontaneous symmetry breaking [f|. 
Obviously, these effects may strongly affect the scatter- 
ing problem. Finally, a challenging problem is to extend 
the analysis to the case of two-dimensional "PT-invariant 
networks. 
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